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Abstract 

In this paper we show that the equihbrium macroscopic entropy of a generic non reversible 
Kawasaki+Glauber dynamics is a non local functional of the density. This implies that equilibrium 
correlations extend to macroscopic distances. 
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I. INTRODUCTION 



In the last two years considerable progress has been made in understanding stationary 
non equilibrium states (SNS). An important result, first established by Derrida, Lebowitz 
and Speer P, and then rederived in a simpler way from general principles in is the 
non local character of the entropy in a boundary driven (i.e. interacting with reservoirs 
having different chemical potentials) simple exclusion system. The question naturally arises 
whether this is a general property of SNS. A different model, the boundary driven zero 
range, does not exhibit this property but one may suspect that this is related to the 
special character of the interaction. All these models are conservative in the sense that the 
number of particles is conserved. 

In this note we present an example of an equilibrium state in which the entropy is a 
non local functional of the thermodynamic variables showing in this way that this is not 
an exclusive property of SNS. By an equilibrium state we mean a stationary state of an 
isolated system or of a system in contact with reservoirs characterized by the same chemical 
potential. Our conclusion is that non locality or macroscopic correlations is more likely to be 
a generic feature of non reversible dynamics. Indeed, in the model considered, the generator 
of the dynamics is not self-adjoint with respect to the equilibrium measure. This means that 
detailed balance does not hold for such an equilibrium ensemble. 

The system considered evolves according to the so called Kawasaki+Glauber dynamics 
with generic transition rates for the Glauber process and will be described in the next section. 
Our analysis is based on a general Hamilton- Jacobi equation satisfied by the entropy in 
equilibrium or non equilibrium states and established in 0, |2|- The H-J equation for our 
model has a very complicated structure but for values of the thermodynamic variables near 
the stationary values can be solved by an iterative scheme of calculation. The entropy is 
non local already in the lowest approximation, i.e. at the level of gaussian fluctuations. One 
should mention that there are special choices of the transition probabilities which make the 

n 

entropy of the model local in spite of a microscopic non reversible dynamics pj, however we 
are interested in properties with some degree of stability. A more detailed analysis of our 
example and of the non locality problem will be presented elsewhere. 
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II. THE MODEL 

In this section we follow closely Q]. The model consists of particles on a lattice evolving 
according to two basic dynamical processes: 

• i) a particle can move to a neighbouring site if this is empty 

• ii) a particle can disappear or be created in a site, the rate depending on the nearby 
configuration. 

The first process is conservative while the second is not. 

Mathematically we consider a family of Markov processes whose state space is X^v = 
{0, 1}^^, where is an integer and Zn denotes the set of integers modulo A^. We shall 
denote with rj a point in the state space, that is a configuration of the system. This is 
therefore given by a function T]{i) defined on each site and taking the values or 1. For each 
N the dynamics is defined by the action of the Markov generator L^r on functions f{r]) 

LnHv) = ^ E ifiv'''^') - fiv)) + E ci^,v)im - m) (1) 
where the addition in Zm means addition modulo N 



r]{k) j = i (2) 



Vij) = { (3) 
1 - r]{i} j = I 

The rates c(i, rf) depend on the values of //(j) with j within a fixed distance R from the site i. 
They are translation invariant, that is there exists a function c(r^) such that c(i, rf) = c{Tirj) 
where {Tkr]){j) = r]{j — k). Let us consider now the unit interval [0, 1) with periodic condition 
at the boundary and a function 7 defined on [0, 1) and taking values in [0, 1]. Let the 
probability measure on the state space of the system obtained by assigning a Bernoulli 
distribution to each site, taking the product over all sites and defined by 

= 1} = 7(^) (4) 
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The main object of our study is the empirical density /if^: 

1 ^ k 

Let us denote by the distribution law of the trajectories /if'(x) when the initial measure 



/^r(^) = ^ E VNHimx - ^) (5) 



is the product measure ly^ . It is possible to show that Q!^ converges weakly as goes to 
infinity to the measure concentrated on the path p{t, x) that is the unique solution of 

d,p = \dlp + B{p) - Dip) 
^P(0,-) = 7(-) 
with 

B{p) = E^Mv)i^ - Vm) (7) 

Dip) = E^MvHO)) (8) 

where i^p is the Bernoulli product distribution with 7(0;) = p. Typically B{p) and D{p) are 
polynomials in the variable p. 

The equilibrium state corresponds to a density p which is the solution of the equation 
B{p) = D{p) that gives an absolute minimum of the potential V{p) = J^[D{p') — B{p')]dp'. 

The above result is a law of large numbers that shows that the empirical density in the 
limit of large behaves deterministically according to equation ©. We can now ask what 
is the probability that /if^ follows a trajectory different from a solution of (jH)) when N is large 
but not infinite. This probability is exponentially small in and can be estimated using the 
methods of the theory of large deviations introduced for the systems of interest in 161 and 

riQ 

developed in [5|, 17( . The main idea consists in introducing a modified process for which the 
trajectory of interest (fluctuation) is a solution of the corresponding hydrodynamic equation, 
and then comparing the two evolutions. 

For this purpose we consider the Markov process defined by the generator 

+ E^ c{z, r/) [(1 - r^(^))e^(*'^) + r/(^)e-^(*'^)] [f{v^) - f{v)] 

with c, 1]^'^ , rf as previously defined and H{t, x) can be interpreted as an external potential. 
The deterministic equation satisfied by the empirical density is now 

dtp = \dlp - dM-^ - p)d,H) + B{p)e^ - D{p)e-^ ^^^^ 
P(0,-) = 7(0 



Given a function p{x,t) twice differentiable with respect to x and once with respect to t 
this equation determines uniquely the field H. The probability that fi^ in the unperturbed 
system follows a trajectory different from a solution of © can now be expressed in terms of 
H and the polynomials B and D. We introduce the large deviation functional 

npt) = \j','j^dtdxpt{i-pt){d^H^f 

+ /o*« /o' dtdxB{pt){l - e^* + //^e^O (11) 
+ /o*« /o' dtdxD{pt){l - e-^* - 

Let pt be a trajectory in the interval of time [0, to] with initial profile 7; assume that the 
initial configuration of the process r/o is fixed and such that //^ (x) 7(x) as ^ 00. The 
large fiuctuation estimate asserts that 

n?(/^r~/^*'^e[0'^o])^e-^^(''') (12) 

where is the probability distribution of the unperturbed process r]t starting in 770- The 
sign ~ has to be interpreted as asymptotic equality of the logarithms. Equation (fT^ is 
a dynamical generalization of Einstein formula for thermodynamic fiuctuations. From this 
equation one sees that the trajectory that creates a certain profile p{x) with highest prob- 
ability is the one that minimizes I{pt) in the set G of all trajectories that connect the 
equilibrium state p to p{x). 
Let us introduce the quantity 

S{p)= MUpt) (13) 

that we shall call the entropy associated to the profile p. One can show, from the large 
deviation estimate ()12j). that in equilibrium 

PiJ(/z^~p)^e-^^('') (14) 

In this way we recover dynamically the usual Einstein formula. Note that the sign of the 
entropy is the opposite of that used in the physical literature. 
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III. THE HAMILTON-JACOBI EQUATION AND ITS CONSEQUENCES 



The functional I{p) has the form of an action integral associated to the Lagrangian 

'C(P,P) = \jodxpt{l-pt){d^Htf 

+ S^dxB{p,){l-e"^ + H,e^^) (15) 

where p) is determined by ()10j). By Legendre transform we can define the Hamiltonian 

H{p, H) = f dx{]-Hdlp + \{d^Hfp{l -p)- B{p){l - exp H) - D{p){l - exp -H)) (16) 
Jo 2 / 

The entropy S{p) then satisfies the Hamilton- Jacobi equation 

^(P,^) = (17) 

This is a very complicated functional derivative equation which however can be solved by 
successive approximations using as an expansion parameter p — p where p is a solution of 
B[p) = D{p) that is a stationary solution of hydrodynamics. For p = p we have g = 0. We 
are looking for an approximate solution of (fTTj) of the form 

^(P) = o / dx dy{p{x) - p)k{x,y){p{y) - p) + o{p- pf (18) 



2 Jo Jo 

The^kernel k{x, y) is the inverse of the density correlation function. It is more convenient, 
as in [2j , to work directly with the correlation function using the Legendre transform of the 
entropy G{h), which is the pressure corresponding to the chemical potential profile h 

G{h) = snp{{h,p)~Sip)} (19) 
p 

G(h) satisfies the dual Hamilton- Jacobi equation 

«(f.ft)^0 (20) 
We are looking for an approximate solution of (|20|) of the form 



G{h) = I dxh(x)p + - I dx I dyh{x)c{x,y)h{y) + o{h'^) (21) 
JO 2 Jo Jo 

where c{x,y) is related to k{x,y) by the following relation 

J c{x, y)k{y, z)dy = 5{x - z) (22) 



By inserting ()2ip in ()20j) we obtain the following equation for c{x, y) 



1 1 

-dlc{x, y) - {di - bi)c{x, y - -p{l - p)dl5{x - y) + bQ6{x - y) = (23) 



where 



and 



bo = B{p) = D{p) = do (24) 
It is now easy to find the equation for k{x,y) by combining ()2Hj) with ()22j) 

- p)dlk{x, y) - bok{x, y) - ]^dl5{x - y) + {d^ - b^)5{x - y) = (25) 

If the entropy is a local functional of the density, k{x, y) has the form k{x, y) = f{p)6{x — y) 
which inserted in fl25|) gives 

m = m-pr' (26) 

and 

feo[p(l-p)]-^-(rfi-foi) = 0. (27) 

Therefore if bo,bi, di do not satisfy this equation the entropy cannot be a local functional of 
the density. 

The equation for c{x, y) simplifies by writing 

cix,y) = p{l - p)6{x - y) + b{x,y) (28) 

obtaining the following equation for b{x, y) 

- ldlb{x, y) + (di - b,)b{x, y) = (p(l - p){b, - d,) + bo)S{x - y) (29) 

Remark that rfi — 61, being the second derivative of the potential calculated in a minimum, is 
positive and the solution of equation is exponentially decreasing. The macroscopic cor- 
relations in the equilibrium states of our system are therefore of short range when compared 
to the stationary non equilibrium correlations of the simple exclusion process considered in 

a. 

In the gaussian process describing, in the same model, central limit type fluctuations 
was studied. When the t ^ 00 limit is taken one finds that the stationary correlations 
of such a process agree, as one could expect, with the macroscopic correlations we have 
calculated in the gaussian approximation. 



IV. CONCLUSIONS 



The calculation presented in this paper supports the conjecture that macroscopic corre- 
lations are a generic feature of equilibrium states of non reversible lattice gases. It shows 
in addition that the Hamilton- Jacobi equation provides an effective approach to the study 
of the macroscopic entropy in cases which are more complex than the conservative lattice 
gases considered in 1^. 
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